Decay of quantised vorticity by sound emission 
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It is thought that in a quantum fluid sound generation is the ultimate sink of 
turbulent kinetic energy in the absence of any other dissipation mechanism 
near absolute zero. We show that a suitably trapped Bose-Einstein con- 
densate provides a model system to study the sound emitted by accelerating 
vortices in a controlled way. 

PACS numbers: 03.75. Lm, 47.32.Cc, 67.40.Vs. 



1. BACKGROUND 

If quantised vorticity is created in helium II at temperatures so small 
that the normal fluid fraction (hence viscosity and friction) is negligible, it 
is found that it rapidly decays* . It is thought that the mechanism which is 
responsible for kinetic energy dissipation is the creation of sound waves, or 
phonond^. It was Nore et aZP who first found that during the evolution 
of a vortex tangle the total kinetic energy decreases and the total sound 
energy increases. By doing a more detailed analysis we revealed that a short 
sound pulse is emitted at each vortex reconnectioiPl The importance of vor- 
tex reconnections was highlighted when it was discovered thatP the vortex 
cusp which is created at each reconnection event relaxes and triggers large 
amplitude Kelvin waves (helical displacements of the vortex core); these 
waves interact non-linearly and generate more waves of higher and higher 
wavenumber. This Kelvin wave cascade 1 ^, which is analogous to the Kol- 
mogorov cascade of classical turbulence, thus shifts the kinetic energy to 
wavenumbers which are large enough that a vortex line can efficiently ra- 
diate souncP. Recently we demonstrated^ the simultaneous occurrence of 
both processes (reconnection pulses and sound radiation). 
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Fig. 1. (a) Radial density distribution n (solid line), including an off-centre 
vortex, for a BEC confined by a tight Gaussian dimple in a weaker harmonic 
trap, (b) Renormalised fluid density (actual density minus inhomogeneous 
background density) for Vq < reveals the emission of dipolar sound waves 
from the accelerating vortex, with amplitude ~ 1% no- 

Unfortunately the study of quantised vorticity in helium II suffers from 
a lack of direct visualisation. Better detection techniques exist for trapped 
weakly-interacting atomic Bose-Einstein condensates (BECs). The key in- 
gredients of the problem (quantised vorticity, sound waves, reconnections) 
are present in both systems. The aim of this paper is to show that a trapped 
Bose-Einstein condensate provides a model system to study sound radiation 
by vortices in a controlled way. Sound radiation in a BEC is thus not only 
interesting per se, but it can give insight into the more difficult problem of 
superfluid turbulence. 



2. MODEL 

In a recent papei^ we have suggested that sound emission by a quan- 
tised vortex can be studied in a controlled way by letting a vortex pre- 
cess within a dimple embedded in a weaker harmonic potential which con- 
fines a quasi-two-dimensional (quasi-2D) atomic condensate, as illustrated 
in Fig. 1(a). If the dimple depth Vq is less than the chemical potential \i the 
sound (which has an energy of the order of [£) escapes, otherwise it remains 
in the region near the vortex and can be reabsorbed. This configuration thus 
allows us to control sound radiation in a sensitive way. 

Our analysis is based on numerical simulation of the Gross-Pitaevskii 
(GP) equation in the quasi-2D limit Fig.^a) shows the density profile of 
the trapped condensate in the presence of a vortex which is initially located 
near the axis of the trap. The vortex precesses around the axis, due to the 
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Fig. 2. (a) Density plot of a corotating vortex-vortex pair in a homoge- 
neous system. Quadrupolar sound emission, of amplitude ~ 0.2% no, is 
generated from the corotating vortex cores (illustrated by plus signs). (b) 
Acceleration of the pair as a function of the radius r (bottom axis), calcu- 
lated from the time-dependent GP equation (black points with error bars), 
a time-independent numerical approach in the rotating frame (solid black 
line), and the analytic prediction^! corresponding to a = K 2 /r 3 (dashed 
line). The acceleration of a single vortex precessing in a harmonic trap 
u> = \f2 x W~ l (/j,/H) is also shown (grey line), as a function of distance from 
the condensate edge (top axis) i?(£) = Rc — r (£)> where the condensate 
radius Rc = 10£. 

Magnus force: the density gradient due to the trapping potential causes a 
buoyancy force which induces tangential motion. What interests us here is 
the fact that the acceleration of the vortex induces emission of sound. 



3. RESULTS 

If Vo 3> n (deep dimple) the vortex moves around the trap in a closed 
orbit, maintaining a mean distance from the axis. The sound emitted by 
the accelerating vortex re-interacts with it and is absorbed. The vortex 
energy does not decay but executes small oscillations which result from the 
beating of the vortex precession and the modes of the trapped condensate. 
If Vq -C u (shallow dimple), the sound radiated by the accelerating vortex 
leaves the dimple, the vortex energy slowly decreases, and the vortex spirals 
outwards to lower densities. The sound waves are emitted in the direction 
perpendicular to the instantaneous direction of motion in the form of a 
dipolar pattern. The precessional motion of the vortex converts the dipolar 
emission into a spiral wave pattern- see Fig. EJb). For a quasi-2D BEC, 
we fincP that the power radiated by the vortex is P = (3mN(a 2 /to), where 



C.F. Barenghi et al. 




t &c) d © 



Fig. 3. Sound burst produced by the close approach of three vortices, (a) 
The trajectory of the vortex-antivortex pair (solid line), initially located at 
(0£, 20£) and (5£, 20£), is deflected as it approaches a single vortex located 
at the origin. The density plot shows the final density distribution of the 
vortices (corresponding to dark spots), (b) The final density distribution, 
shown on a different density and length scale, shows a burst of sound which 
propagates radially outwards, (c) Acceleration experienced by the vortex 
in the pair nearest to the single vortex, for a vortex pair initially located 
at (d,20£) and (d + 5£,20£), where d/£ = (solid line), 1 (dashed line), 2 
(dotted line) and 4 (dot-dashed line), (d) Final radius of the pair following 
the interaction (rescaled by the initial radius of 5£) , as a function of distance 
d. 



Decay of quantised vorticity 



a is the vortex acceleration, u is the vortex angular velocity (induced by 
the trap), N is the number of atoms in the BEC, and the dimensionless 
parameter [3 ~ 6.3 is determined by fitting. Our result is in fair agreement 
with calculations, based on the analogy to (2+l)D electrodynamics^ (where 
vortices and phonons map onto charges and photons, respectively^) and 
classical acoustical, which yield (3 = vr 2 /2, assuming perfect circular motion 
for a point vortex in a homogeneous system. 

In Helium II, the relevant parameter is the decay of vortex line length 
L, typically given by^dL/dt = AirP /[Nmn 2 ln(6/£)], where k = h/m, and 
b is the average intervortex spacing in the tangle. Assuming a generalisa- 
tion of our single vortex power radiation formula to a system of vortices 
with average separation L -1 / 2 , we can cast our result into Vinen's form^' 
dL/dt = — a(K/2ir)L 2 , where a = /3/[vr ln(£L 1//2 )]. This generalisation of a 
numerically-computed expression based on the GP equation yields the cor- 
rect power law decay of vortex line length. Using typical Helium II numbers, 
we obtain a ~ O(0. 1) for the scenario considered here corresponding to a 
vortex decay driven by density inhomogeneity. This is larger than the value 
obtained for the 'free' decay arising from Kelvin wave excitations^, but still 
smaller than the experimentally obtained value for Helium II 9 . To establish 
a link between single vortices driven by trap inhomogeneity and superfluid 
turbulence we first consider the simple example of a vortex-vortex pair. As 
the vortices corotate about their central point, they radiate quadrupolar 
sound waves, which form a double-armed spiral wave pattern, as shown in 
figure Fig.^Ka). Fig. EKb) shows the acceleration of a vortex- vortex pair in a 
homogeneous system and a single vortex in a harmonic trap. Although the 
acceleration is similar in both cases, the sound emission from the pair is less 
due to the quadrupolar character of the emission. This suggests that our 
value for a is an over-estimate for the case of a vortex tangle. The accelera- 
tion of the vortex- vortex pair tails off as the vortex cores merge (r — ► 0), in 
contrast to the analytic prediction of PismerPI. Note that the acceleration 
of a vortex in a harmonic trap is strongly dependent on the trap frequency. 

As an example of a three-vortex interaction, we consider the close ap- 
proach of a vortex-antivortex pair towards a single vortex, as illustrated in 
Fig. Of a). The single vortex and its nearest neighbour in the pair have the 
same polarity. During this interaction, the trajectory of the pair becomes 
deflected (black lines in Fig. 12(a)), while the single vortex makes small devi- 
ations about the origin (not shown). The sharp acceleration of the vortices 
(indicated in Fig. Etc)) induces a sound burst, which propagates radially 
outwards. This is visible in the final density distribution shown in Fig. |3] 
(b). This energy loss results in a reduction in the size of the vortex pair, 
shown in FigEKd). The magnitude of the loss decreases as the initial posi- 
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tion of the pair is moved further away from the single vortex, in line with 
the decreasing acceleration experienced by the vortices (FigJ5Jc)). Note that 
the magnitude of the acceleration induced by the close approach is similar 
to the values studied in the dimple trap. Finally, if the single vortex and its 
nearest neighbour in the pair have opposing polarity, the interaction tends to 
involve a reconnection whereby these vortices form a pair, and leave behind 
the other vortex. 

In summary, we have shown how sound emission can be controlled and 
quantified in experiments on dilute Bose-Einstein condensates. We discuss 
how the information gained can be mapped onto the understanding of the 
decay of superfluid turbulence in the limit of low temperature. However, 
further work is needed to exactly quantify the sound emission for the case 
of many vortices. 
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